Dirac Kirchhoff diffraction theory 
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Kirchhoff 's scalar diffraction theory is applied throughout photon and electron optics. It is based 
on the stationary electromagnetic or Schrodinger wave equation, and is useful in describing inter- 
ference phenomena for both light and matter waves. Here, Kirchhoff's diffraction theory is derived 
from the relativistic Dirac equation, thus reformulated to work on Dirac spinors. The parallels with 
the "classic" scalar theory are highlighted, and a basic interpretation of the result obtained for the 
Fraunhofer diffraction limit is given. The goal of this paper is to emphasize the similarity between 
the two equations on the most fundamental level. 

PACS numbers: 42.25.Fx,03.30.+p, 03.65.Pm 



I. INTRODUCTION 

Scalar diffraction theory is what lies at the basis of 
Fourier optics, and is applied in many different fields in 
optics |i|-|3[ . The formal mathematical theory of diffrac- 
tion was developed by Gustav Robert Kirchhoff in the 
late 1800's. He succeeded in removing the arbitrariness in 
the Huygens-Fresnel theory. Later, Rayleigh and Somer- 
feld removed the mathematical inconsistency formed by 
Kirchhoff's boundary conditions. This led to two distinct 
possible solutions, depending on what boundary condi- 
tions are chosen for the physical model. Kirchhoff's orig- 
inal result was found to be the arithmetic average of the 
two Rayleigh-Somerfeld solutions. For simplicity, we con- 
sider only Kirchhoff's original premise here. 
While investigating electron diffraction through an aper- 
ture, which recently led to the experimental observa- 
tion of electron vortex beams 0-01, the question arose 
whether the simple description of the diffraction process 
by a single Fourier transform was also valid for relativis- 
tic fermions. This goes along with the question of what 
happens to the particle's spin. Both issues can be an- 
swered by solving this problem using the Dirac equation. 
There is work being done on developing the theory be- 
hind electron vortices 0, Q i and there have been previous 
developments with respect to the holographic generation 
of vortex beams Q . A fully relativistic description of the 
creation of electron vortices is useful to ascertain the ex- 
act form of the incoming waves in scattering experiments 
performed using a wave created by a special aperture. 
Usually, scattering theory uses plane waves to describe 
the interaction process, but a vortex cannot possibly be 
simply approximated to a plane wave. A proper form for 
the wave created by binary holograms is required if one 
is ever to develop a fully relativistic description of the 
electron vortex scattering process as it happens in cur- 
rent experiments. 

Explicit calculation of the optical diffraction theorems 



based on the Dirac equation led to a pleasing result: the 
Dirac equation lends itself to a generalization of the clas- 
sical result, containing variables like spin through the in- 
troduction of 4-component Dirac spinors. The formulas 
obtained here can be used to derive more precise calcula- 
tions based on the Dirac equation for a variety of fermion 
wave optics experiments, only by relaxing some approxi- 
mations assumed here for clarity. The derivations in this 
work are based on Hillion's work on the spinor descrip- 
tion of polarized light 10 13], which at its core, uses a 
2-component massless Dirac equation. 
Below, a short derivation of the scalar diffraction theory 
is given containing the essential elements required to gen- 
eralize it for the Dirac case. Afterwards, the exact rela- 
tion between the Kirchhoff and Rayleigh-Somerfeld theo- 
ries will be discussed as it pertains to the Dirac equation. 
Natural units h = c = 1 are used throughout, shorthand 
partial derivative notation is used: dr,., = and Ein- 
stein summation is implied over three (Latin letters) or 
four (Greek letters) dimensions where appropriate. For 
a more rudimentary derivation of the scalar theory, re- 
fer to your favorite book on optics or Hillion's paper on 
the spinor Helmholtz equation [llj , where a 2-component 
spinor form of EM waves is used. 



II. SCALAR DIFFRACTION THEORY 
A. Helmholtz-Kirchhoff integration theorem 

The time-independent Helmholtz equation can be de- 
rived from the stationary Schrodinger equation as well as 
the optical wave equation for the scalar electric field E. 
In what follows, the quantum mechanical notation of ^E* is 
used, to ease the transition to Dirac spinors. The static 
wave equation for an eigenstate ^ with energy E = 
is given by: 



V^*(r) -1- /c2*(r) =0. 



(1) 
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This is the static Helmholtz equation. Setting the right- 
hand side to —S{r — r'), and keeping the part of the solu- 
tion that satisfies Sommerfeld's radiation condition UM, 
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gives following Green's function: 

^ik\r — r' \ 



G(r,r') = 



A-kW 



(2) 



By using the Gauss-Ostrogradski theorem, one can derive 
Green's theorem [15} . 

V 

= §i-en) ■ {U1VU2 - U2VUi)<Prs, (3) 



Where Ui and U2 are two scalar functions that will be 
filled in shortly, is a closed volume with bordering sur- 
face S, and vectors lie on that surface. Filling in 
Ui = -^{r) (a solution of ^) and U2 = G{r,r') into 
([3]), one can obtain the Helmholtz-Kirchhoff integration 
theorem: 



*(rs)e„ • V- 



ik\rs-r'\ 



\rs 



^ik\rs~r'\ 



\rs 



— e„ • V^'(rs) 



<frs. (4) 



This equation gives an exact solution to the wave equa- 
tion ([Ij in any point r' of V, given the value of ^'(rg) on 
the bordering surface S. Being a fundamental result of 
a physically implied Helmholtz equation and the mathe- 
matically rigorous Gauss-Ostrogradski theorem, it finds 
applications in many fields of physics (albeit in a more 
general form). The above can be seen as the mathe- 
matical basis of cosmologist's holographic principle p^ : 
information encoded on a (hyper) surface completely de- 
fines the internal state inside the entire volume enclosed 
by that surface. 



B. Boundary conditions, approximations and 
Fraunhofer diffraction 



The mathematical model of the diffraction setup is 
shown in fig. [T] A point source is located at Pq at dis- 
tance To from an aperture Ai. At point Q G Ai, light 
can be transmitted into the right bounded volume which 
contains a screen at P at a distance s from the aperture. 
We will now calculate the disturbance at point P caused 
by the point source Pq. ^3 is a surface that mathemat- 
ically goes to infinity. The boundary conditions used by 
Kirchhoff are: 

• ^'(r) = dr'^{r) =0 on A2 and ^3, 

• v[/(r) — on Ai, i.e. for r = rg, which repre- 
sents a field originating in a point source Pq which 
illuminates the aperture. 



: Ai s~ 




FIG. 1. Schematic difTraction setup. A point source in Pq 
emits a field which illuminates an aperture A2 at distance ro. 
The wave passes through the aperture [z = 0) at Q and is 
detected on a screen at point P. The outgoing wave makes 
an angle x with the incoming wave. Ai is the unilluminated 
surrounding surface, and A^ is the surface at infinity, s is the 
distance between the aperture and focus point P, where the 
diffraction pattern is calculated. 



Note it is the first of above conditions that breaks down 
mathematically; for an analytic function that has both 
its value and its derivate on a finite region, it must 
be zero everywhere, which cannot be physically true. 
Rayleigh-Sommerfeld diffraction theory sets only one to 
zero, resulting in two solutions that are mathematically 
sound. The problem is choosing the correct boundary 
conditions. In reality, the Kirchhoff approximation pro- 
duces very useful results nonetheless, to a high degree of 
accuracy in the approximations assumed here. 
Some additional approximations are necessary to get a 
useful result: we let fcrg 3> 1 and fcs ^ 1, which means 
that both the source and focal plane are far away from the 
aperture, reducing wavefront curvature effects which take 
more care to be included. This leads to the Kirchhoff- 
Fresnel diffraction formula: 

47r rps 2 

This is valid directly behind the aperture. The quan- 
tity "°'^^"+'^°"^ is called the obliquity factor. 6*0 ~ is 
valid for a point source far away from the aperture, and 
K, Q amounts to a small angle far field approximation. 
The far field approximations and their near field improve- 
ments are discussed extensively in optics literature [it}. 
and not of primary interest here. For completeness, the 
Fraunhofer diffraction formula that one can derive by ap- 
plying these approximations to ([SJ can be written as fol- 
lows: 



^{P) cx JJ< 



--i{k^X+kyY) 



^Ap{X,Y)dXdY, (6) 



with P = {kj;,ky) a point in the focal plane (given in 
Fourier space, which can simply be rescaled to represent 
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real coordinates in that plane) , (X, Y) are the coordi- 
nates on the aperture, and ^'Ap is the wave function in 
the aperture. This formula states that one can easily 
calculate the far-field image (including not only density 
but also phase) if one knows the wave function in the 
aperture. Also note that this theory applies equally well 
to photon waves (described by the wave equation which 
is fully relativistically applicable) and the Schrodinger 
equation (which applies to non-relativistic matter waves) . 
In what follows, the derivation done by P. Hillion (lol[T]| 
in light of the 2-component spinor electromagnetic field 
is expanded to the full 4-component spinor Dirac theory 
for massive particles with spin. 



III. SPINOR DIFFRACTION THEORY 



can construct a solution to the Dirac equation ^ from 
scalar solutions of the Klein-Gordon equation i/;^: 



* = (7^a, - ^k) I 



(12) 



This is a well-known fact. As an extension, which is also 
used by Hillion 11], one can construct a Green's matrix 
solution to 

- d{r - r') = g(r, r') {-i-f^dj + k) , (13) 
from the scalar Green's function ^ as follows: 

g{r,r')=.{ij^dj+k)G{r,r'), (14) 
which is a useful form used in the next section. 



A. Dirac equation 

The Dirac equation for spin-i particles is: 

(i-ff^d^ - m) 4- = 0, (7) 

which is relativistically covariant, and in most common 
usage a 4-component matrix equation for the Dirac spinor 
^P. Applying the complex conjugate Dirac operator to 
the Dirac equation gives an interesting and well-known 
result: 

i-i-r^df, - m) (i7^a^ - m) = 
^ {v'^'d^d^ + m^)^^ 0, (8) 

where i]'^'^ is the Minkowski metric of spacetime. The 
above is nothing less than the scalar Klein-Gordon equa- 
tion for every component of the Dirac spinor indepen- 
dently. This property will be used in the following sec- 
tion. In a stationary situation, one can write ^{r,t) oc 
^'(r)e~*^*, which results in 



{d^dj + ni^) 5- = 0. 



(9) 



One can now define a new variable, k^ = m? -\- , which 
leads to a familiar form: 



{d^dj +k'^)-^ = 0. 



(10) 



The stationary Dirac equation can now be rewritten us- 
ing this new variable as 



(7^9^ +ik)^ ^ 0, 
without loss of generality. 



(11) 



B. Green's matrix 

Using the fact that the solutions of the Dirac equation 
form a subset of those of the Klein-Gordon equation, one 



C. Kirchhoff diffraction 

In this section, a 4-component spinor analogue of the 
Helmholtz-Kirchhoff integration theorem (|4]) is derived 
from the Dirac equation. Multiplying ([7]) with g{r,r') 
on the left and subtracting that from (fT5)l multiplied on 
the right by ^'(r), one has: 



idj [g{r, r')7^*(r)) = ^'{r)S{r - r'). 



(15) 



By using the correct mathematical formulation, one can 
easily prove that the Gauss-Ostrogradski theorem also 
applies to spinor objects [ll|, [11] , in which case one can 
integrate (jlSp over a volume V around the point r' and 
write: 



^'(r') = (0g(rs,r')7"1'(rs)c?^ 



(16) 



where the integration is over the closed surface S around 
V, containing all the vectors rs- 7" is the compound 
7 matrix perpendicular to the surface. Obviously a 7 
matrix cannot be perpendicular to a surface, and this is 
a purely notational convenience. An example of 7" can 
easily be given in the case of a spherical surface, namely 
7'' = sin 9 cos (p'j^ + sin 6 sin ^7^ -I- cos ^7'^, which is noth- 
ing more that the vector transformation to spherical co- 
ordinates. It is not surprising the 7 matrices transform 
as vectors, as they are the basis vectors of the Dirac al- 
gebra. This makes 7" more than only "notational conve- 
nience" , and a mathematical 7 object can easily be seen 
as orthogonal to a surface. From ([2]) one can prove 



-f^djG{r, r')^ik(l + — ] -/^G{r, r')d, 



(17) 



This expression can be used to rewrite (jl6p into a form 
comparable to (HJ: 



■ik nn pik\rs-r'\ 



\rs - r' 



1 



7"*(rs)dVs. (18) 
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The next steps are exactly the same as in the scalar 
theory: Kirchhoff 's boundary conditions and approxima- 
tions are applied, along with: 

• 7'' 7'^: when the beam is directed along the z 
axis, 7*" coincides with 7^^, 



— 1: valid for calculations almost directly 
behind the aperture. 

gives the spinor Kirchhoff- Fresnel diffraction formula 
(compare with 



ir') . JJ. 



Jk(ro + s) 



ros 



(1 + 7^) 7^*(r-s)d5. 



(19) 



D. Fraunhofer approximation 

One can again perform the exact same approximations 
as in the scalar case to reduce p^ . leading to a similar 
form as in ([6]): 

«'(P) oc JJe-^('=--^+'=«^) (1 + 73) -f^^^^J^pdXdY, (20) 

which gives the Dirac spinor form for the far-field image 
of a diffracted z-oriented beam. 



IV. DISCUSSION 

It should not be surprising that scalar diffraction the- 
ory extends cleanly to Dirac 's description of fermions. 
Besides the basic properties of Dirac spinor solutions 
highlighted in Section IIII Al the 2D Pauli algebra of real 
space, generated by {cr*}, is a subalgebra of the 4D Dirac 
spacetime algebra generated by {7''} (both of which obey 
the same anticommutator relationships). The most im- 
portant properties of this result will be discussed: the 
non-relativistic limit and what happens to spin under 
diffraction. 



A. Spin and the non-relativistic limit 

The advantage of working in the Dirac representation 
is the simplicity of its non-relativistic limit. Taking the 
limit — > 0, it is readily known only the top two 
components of the Dirac spinor contribute, and the 7- 
factor effectively reduces to a negligible scaling operation 
for a spin eigenstate. By this, the formula (pO)) reduces 
simply to the classical form (|6]) in the non-relativistic 
case. 

The theory derived here highlights the simple fact that 
for diffraction through a scalar aperture of a paraxial rel- 
ativistic fermion beam, spin is conserved. This is made 
clear by the 7 factor which is multiplied with the incom- 
ing spinor, mixing only the large and small components 



of the same spin elements, the explicit form for the z- 
direction given below for the Dirac representation: 



(21) 



For the other directions, the nonzero elements of the 7- 
factor align with the nonzero elements of the plane wave 
spinor with its spin aligned in that direction. For exam- 
ple, the spinor part of a plane wave eigenstate of = -1-1 
is given by 




/ 1 \ 





(22) 



E+m 

V / 



which, after multiplication with (1 + 7'^)7'^, is still in the 
same eigenstate of E^. The Fourier transform will not 
affect the spinor structure. Hence, the far field spinor 
given by ((20)) is also unmodified with respect to spin af- 
ter diffraction. This result implies that a pure diffrac- 
tion experiment with a scalar aperture will not be influ- 
enced by electron spin, at least in the paraxial region. 
Some mixing may occur far from the beam axis, where 
(1 -|- 7*^) 7'' is not well approximated by (for the z direc- 
tion) (1 + 7'^) 7"^. The generality of the results derived 
here allow to extend the simple approximation (1201) to 
more complicated scenarios, including any type of spin 
polarization. 



B. Kirchhoff and Rayleigh-Sommerfeld boundary 
conditions 



It is well known that Kirchhoff 's boundary conditions 
lead to physical solutions containing a grave mathemat- 
ical inconsistency. By putting both the wave function 
and its derivative to zero on a finite region of space, one 
has inevitably that 4*, an analytic function, must be zero 
everywhere inside the bounding surface (this includes all 
points P) . Rayleigh and Sommerfeld solved this by deriv- 
ing two possible integration formulas, each for different 
boundary conditions: one has only the wavefunction zero, 
the other has only its derivative zero on the "dark areas" . 
Table |T] attempts to give an overview of the conceptual 
differences between the three integral formalisms. The 
Kirchhoff-Helmholtz integration theorems for these cases 
under the present approximations differ slightly: 



Kirchoff: ^ (g {{dr-^)G - ^^^G) dS 
s 

R-S 1: - ^^-^drGdS 

s 

R-S 2: -!- ^{dr-^)GdS 
27r JJ 



(23a) 
(23b) 
(23c) 
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Kirchhoff 


R-S 1 


R-S 2 


d-r'fir) = *(r) = 


9.*(r) = 


*(r) = 








(cos 00 + cos 0) /2 


COS S 


cos 9o 



TABLE I. Comparison between essential differences between 
scalar Kirchhoff and Rayleigh-Somerfeld theories. The first 
row shows the boundary conditions applied on Ai,s in Figure 
[T] The second row shows the Green's function used. The 
third row show the resulting obliquity factor. 



In the Dirac case, the Green's functions become the 
Green's matrices derived from the scalar Green's func- 
tions, through (ITi|) . Note that one can also apply the 
conjugate Dirac operator in that equation, and the re- 
sulting Green's matrix will be antisymmetric to the one 
used here (see [HI), which is expected from the anticom- 
mutator relationship fuUfiUed by the Dirac spinors. 



V. CONCLUSION 

The formal scalar Kirchhoff diffraction theory is shown 
to extend mathematically to relativistic Dirac spinors, 
starting from the Dirac equation. The entire derivation is 



shown, and the result is consistent with the scalar theory. 
The mathematical fallacy present in the Kirchhoff deriva- 
tion is explained, and the effect on the results derived 
here is elucidated. The non-relativistic limit is found to 
reduce to the classical results. The effect of spin is ana- 
lyzed by using the spinor Fraunhofer diffraction equation, 
where it is evident that straight propagation through an 
aperture does not alter the spin eigenstate. More oblique 
incidence and the wave solution at a larger angle from the 
straight line propagation considered here will show dif- 
ferences, governed by the respective directional 7 matrix 
and the incoming spinor wave. 
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